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1. Introduction 

If one uses a lattice regulator to define QCD, renormalisation means that bare parameters, i.e., 
the coupling constant and the quark masses, acquire a dependence on the lattice spacing a such 
that physical observables like hadron masses have a finite continuum limit a — > 0. This should then 
guarantee that all physical quantities have a well-defined continuum limit. 

However, it is not always possible to calculate the observables of interest directly on the lattice. 
Quite often one can only evaluate certain matrix elements of composite operators from which the 
desired physical quantity is derived in a second step. In general, it is then necessary to introduce 
a-dependent renormalisation factors Z for these operators in order to get a finite continuum limit. 
For example, moments of structure functions in deep-inelastic scattering are physical observables, 
but they can only be calculated as products of a (perturbatively calculated) Wilson coefficient and 
a hadronic matrix element of a local composite operator, which is a long-distance quantity to be 
computed on the lattice. 

It is therefore important to investigate the renormalisation of composite operators in lattice 
QCD, and the present contribution describes some of the efforts the QCDSF collaboration has 
undertaken in this field. For the details we refer to Ref. [[I]]. 



2. How to evaluate renormalisation factors on the lattice 

In principle it is possible to calculate the Z factors of multiplicatively renormalisable operators 
by lattice perturbation theory (for a review see Ref. However, perturbation theory on the lattice 
is computationally much more complex than in the continuum and therefore the calculations rarely 
extend beyond one-loop order (see, however, Refs. [^, ||, |5| |6|]). Moreover, lattice perturbation 
theory usually converges rather slowly so that the accuracy of perturbative renormalisation factors 
is limited, even if some improvement scheme is applied. 

Therefore nonperturbative approaches have been developed, in particular methods based on 
the Schrodinger functional (for reviews see Ref. [0]) and the RI-MOM scheme [[§]. It is the latter 
approach in the slightly modified form of the RI' — MOM scheme that was adopted to produce the 
nonperturbative results presented below, where we consider only composite operators constructed 
from two quark fields, Dirac matrices and a few covariant derivatives in between. 

In the Rl' — MOM scheme the basic objects are quark two-point functions with an insertion 
of the operator under consideration at momentum zero. A suitable renormalisation condition is 
imposed, which does not rely on a particular regularisation so that it can be applied on the lattice 
as well as in the continuum. The latter property entails the possibility to use continuum perturba- 
tion theory in order to calculate conversion factors leading from the operators renormalised in the 
RI' — MOM scheme to operators renormalised in, e.g., the MS scheme. 

For a multiplicatively renormalisable operator we express the operator renormalised in the 
scheme ,5? at the renormalisation scale \i as Zj^. e (/i)^(a), where G{a) denotes the bare operator 
on the lattice and the a dependence of the renormalisation factor has been suppressed. For scales 
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satisfying Aq CD <C \i 2 <C I /a 2 the quantity 



(2.1) 

is independent of the scale and the scheme. Here denotes the renormalised coupling in the 

scheme 5? , while j 9 and j8 are the anomalous dimension and the beta function, respectively, 
whose one-loop coefficients are yo and j8o- In the simulations we compute Z^ r ~ MOM and evaluate 

Z RGI = AZ^( M )Z^_ MOM ( At )Z b R l-™( A1 ) . (2.2) 

This expression would be identical to ( [Q| ) if we knew AZ 1 ^ and the conversion factor Z^, _ MOM 
exactly. However, these quantities are calculated in continuum perturbation theory, so they come 
with certain truncation errors and the result will depend on the choice of the intermediate scheme 
J?". What will also matter is the expansion parameter. It turns out to be advantageous not to use 
g m but the coupling g M0M zs as defined in Ref. [§]. 

In bare lattice perturbation theory one-loop results for renormalisation factors are of the form 

ZiL(AOpert = 1 " 7^2 (%M«A0 + A) + 0(g A ) . (2.3) 

As this expansion in the bare lattice coupling g is often poorly convergent, tadpole-improved per- 
turbation theory has been invented, where for an operator with nrj covariant derivatives one has 



Z bare(M)ti 



1 - -fBy ( y Q \n{apL) + A + (n D - \)\% 2 ) + 0(g 4 " 



16?r 2 V 3 



(2.4) 



Here the fourth root of the average plaquette uq = (^trL^n) 1 / 4 is taken from the simulations and 
the expansion parameter is the boosted coupling = g 2 /uq. Combining renormalisation group 
improvement with tadpole improvement one arrives at "tadpole-improved renormalisation-group- 
improved boosted perturbation theory" or TRB perturbation theory for short 

In the above expressions for Z derived from lattice perturbation theory lattice artefacts have 
been neglected and only the logarithmic a dependence has been kept. This is only justified if 
a 2 pL 2 <S 1. Unfortunately, this condition is not always fulfilled in our simulations. On the other 
hand, it is straightforward (though increasingly involved for more complicated operators) to do 
calculations in one-loop lattice perturbation theory with arbitrary values of a 2 pL 2 . One can then 
use the difference between the Zs with and without lattice artefacts to correct the nonperturbative 
simulation results for discretisation errors of 0(g 2 ). 



3. The simulations 



We use gauge field configurations generated by the QCDSF-UKQCD collaborations with two 
degenerate flavours of clover fermions and the Wilson plaquette action for the gauge field. Config- 
urations for 4 values of /3 are at our disposal, /3 = 5.20, 5.25, 5.29, 5.40, corresponding to lattice 
spacings a 0.086, 0.079, 0.075, 0.067 fm. To set the scale for the lattice spacing we have taken 
the value r$ = 0.467 fm for the Sommer parameter, and we use ?"oAms = 0.617 [ pH when the A 
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Figure 1: Z RGI for the tensor current y/ayi V y/ before (left panel) and after (right panel) the perturbative 
subtraction of lattice artefacts. 

parameter of QCD is needed. For each j8 we have between three and five sea quark masses so that 
a chiral extrapolation to vanishing quark mass is possible. 

As mentioned above, we have to compute quark two-point functions with an insertion of the 
operator under consideration, i.e., three-point functions. We restrict ourselves to flavour-nonsinglet 
operators so that only quark-line connected contributions must be evaluated. This is conveniently 



done with the help of momentum sources introduced in Ref. [11]. 



4. Extracting the renormalisation factors 

The simplest procedure for obtaining a value of Z RGI would be to plot the right-hand side of 
Eq. (fO|), i.e., AZ y (jx)!^, _ mqm (ij.)Z^ mom (jj.) versus jj. and to fit a constant to these data in an 



interval of /j. where they form a plateau. Examples of such plots before and after the perturbative 
subtraction of lattice artefacts are shown in Fig. [j] Equivalently one could fit the values obtained 
for Zgjp) = Z^_ MOM ( At )Z b ^ e - MOM ( At ) in the plateau region by AZ y \^- l Z RGl (with Z RGI as 
fit parameter). 

However, there are two effects that jeopardize the reliability of this approach: lattice artefacts, 
which show up at large values of /I and truncation errors of the perturbative expansions in AZ^ and 
^ri'-mom' which become noticeable at small values of ju. They might even conspire to produce a 
fake plateau. Therefore we have tried to incorporate higher terms in the perturbative series treating 
the corresponding coefficients as additional fit parameters. Similarly we have attempted to correct 
for discretisation effects by including a simple ansatz for lattice artefacts. Again, the parameters 
in this ansatz have to be fitted. We fit the data for all four /3 values simultaneously and only the 
quantities Z RGI , our final results, depend on j8, the other parameters do not. Two examples of such 
fits are shown in Fig. ||. For further details concerning the fit procedure we must refer to Ref. [jj]]. 
The results determined by this fit procedure will be called fit results in the following. 

Unfortunately, these fits work only for the perturbatively subtracted data. However, for techni- 
cal reasons, we have applied our subtraction only for operators with at most one covariant deriva- 
tive. Hence we must apply a different procedure for operators with more than one derivative: We 
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Figure 2: Z (perturbatively subtracted) for the tensor current (left panel) and ^7(1 ^4} V (right panel) as 
a function of the renormalisation scale. Also shown are the fit curves used for the determination of Z RGI . 

read off Z RGI at \i 2 = 20GeV 2 and take as the error the maximum of the differences with the re- 
sults at ji 2 = lOGeV 2 and /I 2 = 30GeV 2 . These results will be called interpolation results in the 
following. Obviously, this method can also be applied to subtracted data. The errors assigned by 
the fits appear to be seriously underestimated since they are mainly determined by the statistical 
uncertainties. Therefore it seems to be more reasonable to finally use only the errors from the 
interpolation method as these take into account also some of the systematic effects. 

5. Results 

Due to lack of space we have to restrict the presentation of our results to operators without 
derivatives, i.e. the "currents" 

G s = ud, G p = uy 5 d , (f^ = uy^d , &^ = uy^ysd , ff T ^ = ua^d , (5.1) 

and the quark wave function renormalisation constant Z q . 

In Fig. H we plot the results at /3 = 5.40 extracted from the perturbatively subtracted data, 
both by interpolation and by means of the fit procedure, and the interpolation results based on 
the unsubtracted numbers as well as one-loop perturbative estimates. Ideally the nonperturbative 
results should agree within the errors. In reality, this is not always true. Note, however, that the 
errors of the fit results only account for the (rather small) statistical uncertainties of the raw data 
while the errors of the interpolation results are dominated by systematic effects. The one-loop 
perturbative estimates are larger than the nonperturbative values, but tadpole improvement works. 
TRB perturbation theory, on the other hand, leads to further improvement only in a few cases, for 
some operators it is even worse than ordinary tadpole-improved perturbation theory. 

In Fig. H we compare our fit results for the operators without derivatives with the one-loop and 
two-loop perturbative estimates, again for j3 = 5.40. The numbers from bare lattice perturbation 
theory, represented by circles in the figure, exhibit the expected behavior: The two-loop results 
come closer to the nonperturbative numbers than the one-loop estimates, though only slightly in 
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Figure 3: Results for operators without derivatives at /3 = 5.40. The filled symbols correspond to our fit 
results (circles), interpolation results based on subtracted (squares) and unsubtracted (triangles) data. The 
open symbols represent estimates from bare perturbation theory (circles), tadpole-improved perturbation 
theory (squares) and TRB perturbation theory (triangles) based on one-loop calculations. 




Figure 4: Results for operators without derivatives at j3 = 5.40 compared with one- and two-loop lattice 
perturbation theory. The crosses correspond to our nonperturbative results obtained by fits of the subtracted 
data. The open symbols represent estimates from bare perturbation theory (circles), tadpole-improved per- 
turbation theory (squares) and TRB perturbation theory (triangles) in the one-loop approximation. The 
corresponding estimates based on two-loop calculations are shown by the filled symbols. 
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the case of the tensor current. Except for the tensor current, tadpole improvement works also in 
the two-loop approximation moving the perturbative values closer to the nonperturbative numbers. 
However, the results from TRB perturbation theory, shown by triangles, do not differ much from 
the values found by tadpole improved two-loop perturbation theory. 

6. Concluding remarks 

The RI' — MOM scheme has been established as a method for nonperturbative renormalisa- 
tion that can (relatively) easily be implemented for arbitrary lattice fermions. Momentum sources 
allow us to deal with all operators in a single simulation and at the same time to achieve small 
statistical errors, but the required computer time is proportional to the number of momenta con- 
sidered. One of the largest sources of systematic uncertainties are discretisation effects: Here the 
perturbative subtraction of lattice artefacts has proved very helpful. Continuum perturbation the- 
ory is needed for the conversion to the MS scheme. Obviously, one should use as many loops as 
are available, but additional improvements are possible through the careful choice of intermediate 
schemes and expansion parameters. Still, it seems that the available perturbative results cannot 
describe the scale dependence below the (surprisingly large) value of about 5GeV 2 , as seen, e.g., 
in Fig. [T]. Remarkably enough, there are now renormalisation factors calculated in two-loop lattice 
perturbation theory. Unfortunately, it is difficult to predict their accuracy without comparing with 
nonperturbative results, but (at least in most cases) improvement seems to work. 
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